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The current algebra for gauge theories like QCD at finite temperature and density is studied. We 
start considering, the massless Thirring model at finite temperature and density, finding an explicit 
expression for the current algebra. The central charge only depends on the coupling constant and 
there are not new effects due to temperature and density. From this calculation, we argue how to 
compute the central charge for QCD4 and we argue why the central charge in four dimensions could 
be modified by finite temperature and density. 

I. INTRODUCTION 

Anomalies appear at quantum level because some classical symmetries are not preserved. There are some dramatic 
experimental consequences from the anomalies as, for example in the 7r° — > 2j decay [J. 

Although there are finite temperature and density corrections associated to the ultraviolet sector, the most important 
contribution to those corrections, probably, will come from nonperturbative effects in the infrared sector 00 01 [HIE 0- 
This last fact is particularly relevant when phenomena such as hadronization and deconfinement are considered. 

If, formally, we consider Z as the generating functional 



Z 



VA, 



- slA] Z F [A], 



describing a general gauge field theory, where S[A] is the kinetic term for the gauge fields and Z F [A] is the generating 
functional for the fermionic fields minimally coupled to the A^, then the fermionic vector current in the chiral 
limit is defined as 



6Z F [A] 
6A» 



(1) 



where J M = ip^^ip. 
Then, the relation 



8 2 Z F [A] 



S 2 Z F [A] 



5An(x)6A v (x!) 
J M (x), J v (x ] 



5A u (yL)o~A^(Tt) 



(2) 

provides, in principle, an explicit calculation procedure for the current algebra. 

This last statement needs some clarifications. Indeed, in the past many people computed the current algebra @ 
using canonical methods and point splitting regularization. However, the novelty of our approach is that, firstly, we 
compute the fermionic determinant at finite temperature and density and, later, we compute and show the universal 
character of the current algebra. 

The direct calculation of the current algebra in four-dimensions involves many technicalities and the appearence of 
divergences. It is, therefore, not easy to recognize the temperature and density corrections from such direct calculation. 

These considerations motived us to start considering, initially, a more simple problem in order to get some physical 
insight for the four-dimensional case which will be considered later. 

Thus, we will start considering the massless Thirring model at T = and we compute explicitly the central charge 
From this calculation, we will be able to extract also the current algebra structure for finite temperature and 
density and we will show that the central charges for T = (ct—o) and T^O (cr^o) satisfy 



ct=o = ct^o, 



(3) 



i.e. finite temperature and density do not induce corrections to the central charge. Although the temperature 
independence could be argued given physical arguments, the density independence is nontrivial as we will discuss in 
the conclusions. 
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This letter is organized as follows: in section 2, the current algebra for the Thirring model at T = is considered. 
In section 3, we discuss temperature and density effects in the Thirring model. In section 4 we argue how to extend 
the previous results to higher dimensions. Finally, in section 5 we give the conclusions and outlook. 



II. CENTRAL CHARGE IN THE THIRRING MODEL 



In order to introduce the calculation method, let us start discussing how to compute the current algebra for the 
free case. 

One start considering the lagrangian 

C = ${ip-gA)^, (4) 

where is an auxiliary field which will vanish at the end of the calculation. 
Then, the partition function in the euclidean space is given by 



From JSJ), the vector current becomes 

5Z A 



6A» 



VipVip e~ s[A] . (5) 



= / VipVip tj}{x)^{x) e- s[A] 

A=Q J 

= 9 {M x )) A=o > (6) 



and, therefore, the Ward identity becomes 

9 t ,(J fi (x)) A=0 = 0. (7) 

Although in this free case, this identity is satisfied identically, in the interacting case this conservation law could fail. 
Indeed, the procedure discussed above implies to compute the fermionic determinant and, therefore, a regularization 
scheme must be used. 

If one consider an interacting quantum field theory involving classical gauge and chiral symmetries, then one cannot 
retain simultaneously both symmetries. For example, if we regularize the fermionic determinant with a gauge invariant 
regularization scheme then, the chiral symmetry is explicitly violated and vice versa. 

Some questions, relevants in this context, are for example, which are the implications for the current algebra 
calculation?, will be this algebra modified due to temperature and density effects?, which are the modification in four 
dimensions? and so on. 

From JHJ, one could compute the current algebra by means of 

62Zl ^ =([J,(x),J,(y)]), (8) 



SA^i^SA^iy) 5A v (y)5A»(x) 

using different alternative regularization procedures and not necessarily the splitting point regularization. 
We will consider below the Thirring model 

C T = $ip ^-g 2 (^) , (9) 

and will compute the associated current algebra. 

In order to do that, it is more convenient to use the linearized form of JSJ, v.iz. 

C A = $(ip -gA)^ + - A A\ (10) 

which is equivalent to 10 after integrating A^. 
This last fact implies to change A^ into 

A^x) -> 2g^{x) 1 ^{x) = 2gJ li (x), (11) 

and, as a consequence, (0 and ((Tn|) are equivalent. 
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In this interacting case, the vector current in the sense of © is 



SZ A 



5A^ 



= l(A fl )-g(J fl ), (12) 



where the first term in RHS in (|12|l is the correction to the conservation law Q) due to the interaction. We would 
like to emphasize, however, that the notation (A^) is formal because the expectation value is computed integrating 
the fermionic fields. However, this notation will be useful in the next formula where -by using an appropriated trick- 
the integral in A^ will appear explicitly. 

Although classically the relations (A^) = and 9 M ( J^) — hold, at the quantum level the situation is quite 
different 

Indeed, in order to compute the current algebra one use the identity 



A v {x ) - 2gJ v (x ) 



(a^A^x)) = J V^ViJjVA A^(x)A v (x) e~ s W 
V^Vifi e~ ST I VA [A^x) - 2gJ lx {x)) x 

4 g 2 J V^Vi) J^x)J u {x) e- ST . (13) 

The first term in the RHS is just a gaussian integral as can be checked using the change of variables 

Afj_ — > A'^ = A^ — 2g J^, 

and, therefore, one obtain 

VA' A'^x)A' v (x') e-i / *** A " (1) T = 2S^S 2 (x - x') (1) T , (14) 
where in the last step (1) T is a notation for 

(1) T = / V^Vijj e~ ST . 



The second term in the RHS i|13|) . contains the product of fermionic currents evaluated at two differents points. 
However, as we are interested in the current commutator, the term (1) T is cancelled out and, therefore, the only term 
that finally remains is the commutator < [J M (x), J v (x )] >. Notice that we arc unable to calculate directly this term. 
However, the LHS of 113|) can be computed through the fermionic determinant 

/ VAA,(x)Mx')e-^f d2x l(^ D - lF ' P ) +A2 ]. (15) 

In the momentum space (|15l) becomes 

/ *!L e -«(-*') / VAAj-k)A^k)e'^ d2k ' A ^ k '^ G ^ A ^ = 
J \2^Y J 

= J d 2 kG-l{k)e-^-*'\ (16) 

where the inverse of G^^k) is 

= 2 {-^E - + ^} ' (17) 

Collecting terms, one find that (|13fl becomes 

Ag 2 (J^x)J v (x')) + 2S^S 2 (x - x') (1) = 

^G-l(k)e-^~* '). (18) 
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Thus, once 118fl is obtained, then the current algebra calculation is straightforward . Indeed, from (|18(l one see that 
the current-current commutator is 

([J^x),^!)]) - 

= lim [ ^G-l(k)e-^-* i)sin(fc r). (19) 



In order to compute the RHS in l|19|l . one note that in the limit r — > + becomes 

lim /" rffc , o sin(fcor) = 



2\ 4 



where the property 



= /,,^f ^( V /,77i = -/M (->()! 



lim ifi(z) = ij — , (21) 



has been used. 

The calculation of the diagonal elements is straightforward and it gives 



lim f e?fco ^"^° sin(fc r) = 
r->0+ J k z 

/k\k\ 
dko o sin(fcoT) = 
k z 



lim / dk Q sm(k T) = 0. (22) 



Therefore, the current algebra for the Thirring model becomes 



where the central charge ct is 



([J (x 1 ),J 1 (x' 1 )}) = ^-5'(x l -x' 1 ), (23) 



CT 



1 



(Ml 

2tt 



One should note that the limit g — » is smooth and coincides with the central charge for free fermions. 

However, there is a divergence for ct in g 2 — ^. For this value of g 2 our calculation is no longer valid. Indeed, when 
g 2 = ^ one of the eigenvalues of G^ik) becomes negative and, therefore, the calculation of the gaussian integrals 
becomes meaningless. 

Therefore, the values g 2 > \ will yields necessarily inestabilities and the Thirring model is only properly defined 
for ct < 0. 

III. INCLUDING TEMPERATURE AND DENSITY 

In this section we will explore the effects on the current algebra due to finite temperature and density. 

By convenience we will work in the Minkowski space. In this case the Thirring lagrangian @ is changed by 

Ct-H^PH, (24) 

where [i is the chemical potential and the temperature effects are incorporated in the fermionic modified propagator, 
v.i.z. 

S(k)~> S(k') T=Q + S(k') T ^ , (25) 
where k' v = k u + /i u v and u v being the four velocity with respect to the frame where the heat bath is at rest. 
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In (|25|l S(k)T=o is the standard Dirac propagator for massless fermions, but now including the chemical potential 
contribution, i.e. 



(k + n u) 2 lit' 



while the temperature correction is 

S{k) T ^o = 2iTT{^ + n F {k ■ u) 5[{k + fj, u) 2 }, (27) 
where np(k ■ u) is the Fermi-Dirac distribution given by 

1 



e f)\k-u\ _|_ I ' 



In the fluid rest frame, the vector u is w M = (1, 0), and in such case, the Fermi-Dirac distribution depends only on fco- 
Contrarily to the zero temperature case, the effective action contains contributions to all orders |I2| and, therefore, 
one could think that the finite temperature contributions could induce corrections to the central charge in the current 
algebra. 

In order to explore this, let us consider the effective action up to the second order i.e. 



Vn : / ^A^k)W^{k)A v {-k), (28) 



where the polarization tensor is given by 



g 2 f d\ 



l^S(q)-y v S(q - k) 



(29) 



= n$ (k) + n w (fc) + n$ (fc) + n.l. , (30) 

where N.L. denote the non-local contributions that will be computed in the appendix but, they do not contribute to 
the central charge calculation. 

The equation (|28|l yields three different contributions to the effective action. The first one, formally, is the analogous 

to the T = case including now the chemical potential in the poles of Ilj2 . Naively, the presence of the chemical 

potential will produce a non-trivial result. However this is not true. Indeed, due to Lorentz invariance 11$ can be 
decomposed according to 

= a r7 M „ k 2 + b k^ky. 
Additionally, the gauge invariance should imply that 

fc"n$ = 0, (3i) 

and then from H3I|) a = —b. 

The explict form of a is obtained by computing fl^ ^ i.e. 

a = n. (32) 

A straightforward calculation for LT^ ^ = II yields 

9 2 
fl= 27' 

and, therefore, the chemical potential does not play any role in the first term of the effective action. 

This last fact is, of course, a two-dimensional accident. Physically, this means that there is no compact matter in 
two- dimensions. 

The remaining diagrams could give temperature and density contributions. The sum of the diagrams with one 
thermal insertion (fig. 1) are 



k ft k k 

Q-k 1-k 



FIG. 1: Diagrams with one thermal insertion 



n«(*) = ng?>(*) + n<»>(fc) 

cPq 



2 

t. 

2 



(2n) 2 
d 2 q 
(2n) 2 



Tr 
Tr 



In S(q) T =o 1v S{q- k) T ^o 
Iti S(q) T ^o 7" S(q- k) T =o 



and explicitly each diagram is 



nig 



d 2 q 
{2nf 



Tr 



(q + (i uY + it 



x n { (/3\u ■ (q - k)\) 6{(q - k + ^ u) 2 ), 



and 



nig 



d 2 q 
(2n) 2 



Tr 



7m + M 



— ty + lA-ft) 



((q — k + fi u) 2 + ie 



x m(0\u-q\) S((q + Liu) 2 ). 



(33) 



(34) 



(35) 



Note that both integrals are related. In fact, since they are convergent we can shift variables. Indeed, the changes 
g->g-fc + /i«, for the first integral in (|34|l and q — > q + ^ u for the second one Ij35(l imply that the previous integrals 
become 



■nig 



d 2 g 



Tr 



K + fl 



and 



Trig 



x n f (/3|u • g - <% 2 ), 



(36) 



Tr 



lii\ 



(37) 



((g - k) 2 + ie 
x u f (/3|u-g-/x|) (5(g 2 ). 

Note now that, the properties of gamma matrices imply that 

Trlj^vp] = Tr[7 M ^^. 

Therefore, the diagram corresponding to (|34(l is exactly the same corresponding to (|35J) by changing k 
Equivalently, the diagram l|36|) is the same as 

133 (n^„ 2) (fc)) changing k -> -k. 
In order to compute one of these diagrams, firstly, one can decompose the tensorial structure in ( I37f) as follows: 

nW(k) = a ( 12 V+/3( 12 )^ +7 (12) %^, (38) 

where the notation At a Bp\ stand for the symmetric product of A and B. One should note that any 2-rank symmetric 
tensor can be written as a linear combination of ry M „, k^k v and k(pU v \, if fc M and are independent. 



7 



It is straightforward to prove that the components a, (3 and 7 can be written in terms of Lorentz invariants as 

Wk v 



a W = n (i2)^ _ n (i2) 



fc 2 ' 
1 



(fc ■ u) 2 — k 2 



{k-u)U^k^u^ + 



,(12) 



2 V(^-u) 2 -fc 2 



fc 2 fc-u 



(39) 



(12) (12) 

Therefore, the only relevant integrals that must be performed are the corresponding to 11^ ' if v , IL^' k^k v and 

(12) (12) 

Straightforwardly one can prove that n^^r^ = 0. The product n^&^A^ is 

2«5 2 J [2(3 ' *)((« - fc) • fc) - fc 2 (9 - fc) • «] x 



,n f (/3(|w g - M |))% 2 



(<7 — fc) 2 + it 
d 2 q 



2{q ' k) k -n f ((3\u-q-n\)5(q 2 ), 



-2nig^ 



k) 2 +ie 
(k-q) ni (P\u-q-fx\)5(q 2 ). 



(40) 



This trace can be recast in a more transparent fashion by rescaling the integration variable (3q — > q. 
This trace becomes 



TlffikWc" = -2mg 



■ 2 k a v a 



(41) 



where 



J^^Qa n f (\u ■ q - f3fi\)5(q 2 ). 



Finally, let us point out that the vector v a has components only along the direction defined by u^. That means 
that 



n^fc" = -4tu^ J ^L(q ■ u) n f (\u ■ q - fa\)6(q 2 ). 
In order to compute the coefficient in front of k^u v \ one proceed as follows: firstly, one compute 



(42) 



d 2 q nf((3\u ■ q — fi\)S(q 2 ) 



(2tt) 2 (q - k) 2 + it 
{(q-k)((q-k)-u) + (q-u)((q-k)-k)- 
(fc • u)((q - fc) • q)}, 



Amg2 J w? iq - u 



2(q -k)-k 2 
(q — k) 2 + ie 



n f ((3\u-q-n\)5(q 2 ) 
d 2 q 



(2ny 



(q ■ u)n { {(3\u ■ q - fi\)5(q ), 



(43) 
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FIG. 2: Diagram with two thermal insertions 



then, this last integral can be recast by introducing the same rescaling of q as before. 



-4-irig 2 f d 2 q 



V J (2tt)2 

The coefficients a^ 12 \ ■ ■ ■ can be evaluated directly, indeed. One find 



(q ■ u) m(\u ■ q - (3fi\)6{q 2 ). 



a («) = 2 7ng 2, ^X(Pfi,u), 



0(12) 

7 (12) 



0. 



where 



Therefore 



d 2 q 



13k 2 

2TTig 2 ^X{(3n. /U ). 



(q-u) n { (\u- q - f3^\)6(q 2 ). 



,k ■ u 



n^(fc) = 27ng 2 'j^-X(Pfi, u) ; ; 



(44) 

(45) 
(46) 
(47) 



(48) 



From this last result we see that this diagram is odd respect to the change k — * —k. Therefore, the sum of diagrams 
and iJSSJl is zero 



n(i 2 )(fc) + nW(fc) = n£?>(*o + ng?>(-*) 
= n^)(*)-n^)(fc) 
= o. 



(49) 



Finally, we must compute the diagram with two thermal s, (fig. 2) given by the expression 

d 2 q 



n$(*) = 2^ 2 /^ 



it turns out that 



x nf[/3\u ■ (k -q)- fi\] n { [(3\u ■ q - n\]6(q 2 ) 

<5(( fc - qf)[iM - k)„ + q v {q - % - v^(q - k ) ■ <?]• 



nS(*) = o 



(50) 



Once these diagramas are computed, then the effective action becomes like (J2S), but with poles in ko shifted to 
fc - /U. 

In the case of gauge theories like QCD, the cancellation of these thermal diagrams is a general and mandatory 
property if the gauge invariance is preserved. 

The calculation of the current algebra for the Thirring model is straightforward and, as a consequence, there are 
no modifications due to finite temperature and density effects in two dimensions. In the next section we will come 
back to this conclusion when we go beyond two dimensions 
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IV. FINITE TEMPERATURE AND DENSITY BEYOND TWO-DIMENSIONS 



In this section we argue how to extend our results to the four dimensional case. One start considering the lagrangean 

C = tp{ifi -gfS -m)i>- -F 2 , (51) 



where the effective gauge field B is 



B„ = A^ + - Vfl , (52) 



the velocity of the fluid and m an ultraviolet regulator. 
The effective action for this case is 

d 4 k 



.9 



S *ff = I 7^Mm^{k)B v {-k), (53) 



where the polarization tensor 

n ptf = nw + n#+n$, (54) 

is now quadratically divergent and and it needs to be renormalized. The contributions in (|54|l include the standard 
self-energy gauge field diagram plus the contributions with one and two thermal insertions. 

In analogy with our previous calculation, one could conjecture that these thermal contributions should vanishes 
due to gauge invariance of the effective action (|53[1 . However, (|38ll does not remain valid in the 4-dimensional because 
we have to introduce a new contribution proportional to it M u M 13]. 

However, the contribution n^°J could contain chemical potential contributions. Actually, this effective action written 
in terms of the original fields becomes 

S eff =g 2 J d 4 kA^k)H^ v A v {-k) 

+ 2g J d 4 k A^k)!!^^ v v +g J d A k v^{k)Ii^^v u . 

(55) 

The renormalized polarization tensor is 

n^-(9 2 ^-9^)n(°)( g 2 ), (56) 

where H^°\q 2 ) is a finite contribution given by 



n<°V) = f 1 dx x(l - x) In 

2t>" Jo 



m 2 — x{l — x)q 2 



. / 4 m 2 arctan( , q Y 

= — log m \-q+ ,,■> ■> 

9 \ \J 4 772 — q 2 



(57) 



where m 2 is a cut-off which is necessary in order to have a well defined integral. 

The diagram implies, as usual, charge renormalization where the effects induced by density and temperature appears 
in the qo component if the substitution qo — fi is understood. 

Thus, one can conclude that the effective action can, of course, depend of the temperature and density and that 
the gauge invariance is a preserved symmetry 

In the case of the current algebra calculation the situation is different. Indeed, corrections to the central charge 
come from the term k^k^ and as this term is multiplied by Ii{q 2 ), then the central charge would be modified by 
density and finite temperature. 

We will analyses this problem elsewhere. 
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APPENDIX A: ONE-LOOP NON-LOCAL CORRECTIONS TO THE EFFECTIVE ACTION 

At the one-loop order, the polarization tensor in components can be written as follows (v.i.z. (1291) 1: 

n°°(fc) = n^(k) = 2 J ^( q ^ k + q r + q \ q + ky)^^^ + (2 m )n F (\k\)S(k^ x 

' - + (2ni)n F (\k + q\)5((q + k) 2 )) (Al) 



{ q + k) 2 +ie 



and, 



Uf(k) = Uf(k)=2 |^(qO (fc + g) i_ g i (g + fc) o ) ^_l_ + (2 „ )nF( | fc | )( 5 (fc2) ^ x 

- + {2m)n F (\k + q\)S((q + fc) 2 ) ) . (A2) 



(q + k) 2 +i 

Our next step is to compute the diagonal thermal corrections of W iV . The first correction comes from 

Am J (2n) 2 (k0) 2 -(k^ + l e Siiq + k) )HF{lk +Q l} - (A3) 



However, before to computing l]A3fl . it is convenient to use the identity 

S((q + kf) = 2|fcl X +gl| (S(q° + fc° - (q 1 + k 1 )) + S(q° + k° + (q 1 + k 1 )) . (A4) 
Therefore, (|A3J| becomes 

2q 1 + (fc 1 =p fc°) 



V f d 2 q V 

(2tt)2 (2gi + (fci 



^ J (2tt) 2 (2q 1 + (fc 1 T&°) + 

xe(fc° + q°)%° - fe° ± (g 1 + fc^Jnjfl? + fc°|). (A5) 
We will make use the Plemelj's decomposition to compute these integrals: 

= V ( J — ie(a)ir8(x — Xq) 



Ct(x — Xq) + it \X — Xq 

where V mean the principal value and the e-prescrition is assumed. 
Using this identity, l|A5(l becomes 



- 2m (j^- + m [S(k° - fc 1 ) + 5(k° + fc 1 )]) J ^-e^ 1 )^ 1 + k^npQk 1 + q^). (A6) 
The next correction to the diagonal terms is calculated in a similar way. Indeed, if we use 

% 2 ) = ^(% + ? 1 ) + %°-<z 1 )), 



One find that 



2m (-ZjL+m [S(k° - fc 1 ) + S(k° + fc 1 )] ) x 



x/^^W+^Vd? 1 !). (A7) 
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Finally, the last contribution including two thermal insertions is 

2(2nif ( 4\i (gV + k°) + <zV + fc 1 )) x 



. (2^) 2 

x WTTO W + k °- & + + + k ° + & + x 

x (S(q° - q 1 ) + S(q° + q 1 )) n F {\q°\)n F (\q° + k°\). (A8) 

After collecting terms, the non-local corresctios to the diagonal terms is given by 

Ilf(k) = 2n 2 (d(k° + k 1 ) + S(k° - k 1 )) (A9) 

^VetfWq 1 + k l )L F {\q l \) + n^q 1 + fc 1 !) - 2n F (\q 1 \)n F (\q 1 + fc 1 !)) (A10) 



(2tt)' 



in full agreement with Das and Da Silva Of course, these contributions do not modify the central charge. 
In the same way, one verify that 

Uf(k) = 2n 2 {S(k° + k 1 ) - S(k° - k 1 )) J AefeW + fcV 

x {n F (\q 1 \)+n F (\q 1 + k 1 \)-2n F (\q 1 \)n F (\q 1 +k 1 \)}. (All) 
All these results are proportional to the integral 

h(k 1 ;P) = j -^WH^ + q 1 )* 

x {npitfl) + npdk 1 +q 1 \)- 2n F (|< ? 1 |)n F (|fc 1 + q^)}} . (A12) 

Although l)A12jl is a cumbersome expression, the limit [3 — > oo can be easily computed. Indeed, since I2 — > when 
f3 — ► 00, then in this limit 

coshf ^HfVl A 

dq 1 e^Hk 1 + q 1 ) . — — ^ . ^, — — v . (A13) 



(2^7^^ cosh^ '^l+'f^ 1 ' ) +co S h(/3 ^'-'f^ 1 ' ) 

If we assume that p 1 > and one consider three regions on the real axis, namely, IZi = (—00, —p l ), IZ2 = {—p, 0) 
and IZ3 — (0, 00). Then in 7?.! and IZ3 the product e(g 1 )e(fc 1 + q 1 ) is positive and in 1Z 2 is negative. 
Thus in these regions we have: 

\q 1 \ = -q\ Ifc 1 + q X \ = -(fc 1 + q 1 ) en^i 
\q l \=-q l , Ifc 1 + q 1 ] = (k 1 + q 1 ) enTZ 2 
\q 1 \=q\ Ik 1 + q 1 ] = (k 1 + q 1 ) en Tl 3 



and, therefore, 



Iff ,1 cosh(/3fcV2) f , cosh(/3(fc 1 +qV2)) 



(2tt) 2 \y WlUTC3 cosh(/3fc 1 /2) + cosh(/3(fc 1 +gV2)) J TC2 cosh(/3(A;9 + g 1 / 2 )) + cosh(/3fc72) 
However, in the last integral one use the identity 

cosh(/3(fc 1 + ^72)) = i cosh(/3fcV2) (Au) 



cosh(/3(fc 1 + q 1 /2)) + cosh(/3fcV2) cosh(/3fcy2) + cosh(/3(fc 1 + gi/2)) 

and, therefore, one obtain that 

2 (2tt) 2 ^ /Sy.^, coBhGS^ + coshfx),' 
where in the last integral the change of variables x = I3{k l + q 1 /2) was performed. 
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Finally, I2 is 



h{k l ;(3) = -7^2^(1 " ootgh(/3fcV2)) = -e 2 1 , (A16) 



where at low temperatures its behavior is given by 



in full agreement with the expected result [l^ . 
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